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Abstract 

Numerical simulations of full QCD on anisotropic lattices provide a con- 
venient way to study QCD thermodynamics with fixed physics scales and 
reduced lattice spacing errors. We report results from calculations with two 
flavors of dynamical staggered fermions, where all bare parameters and the 
renormalized anisotropy are kept constant and the temperature is changed 
in small steps by varying only the number of time slices. Including results 
from zero-temperature scale setting simulations, which determine the Karsch 
coefficients, allows for the calculation of the equation of state at finite tem- 
peratures. 
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I. INTRODUCTION 



Lattice calculations provide a method to study the properties of the quark-gluon plasma 
(QGP), which is considered to be the phase matter existed in at the extreme range of 
temperatures microseconds after the big bang. To understand the basic character of the 
QGP we need to determine the equation of state (EOS) for the system in a regime of 
strong gauge coupling for which a non-pcrturbativc scheme for calculation is most adequate. 
Currently in the lattice formulation of QCD two different approaches to that problem could 
be adopted. One of them is the derivative (operator) method which requires the knowledge 
of the asymmetry coefficients, or Karsch coefficients [1,2]. These coefficients have been 
evaluated only perturbatively [3,4] since their non-pcrturbative values are not trivial to 
calculate in practice [5]. The second method is the integral method [6,7], which does not 
require the values of the Karsch coefficients, but has the disadvantage that for a given quark 
mass at a single temperature a number of different simulations are required, and in addition, 
there exist the problems of scaling violation. In our study we avoid these disadvantages by 
choosing the derivative method implemented for anisotropic lattices in combination with a 
fixed parameter scheme described below. 

The anisotropic formulation of lattice QCD has certain advantages regarding the study 
of the equation of state at various finite temperatures. Finite temperature field theory has 
a natural asymmetry which makes the anisotropic approach useful to reduce the lattice 
spacing errors associated with the transfer matrix at less cost than is required for the full 
continuum limit [8] . Through the introduction of anisotropy on the lattice one can make the 
temporal lattice spacing, at, sufficiently small so that by varying the number of time slices, 
N-t, the temperature can be changed in small discrete steps. 

To study the thermodynamic properties of the quark-gluon system we simulate QCD with 
two fiavors of dynamical staggered fermions on anisotropic lattice. The fixed parameter 
scheme we employ to avoid scaling violations is the following. All the bare parameters 
of the simulation are kept constant and only the temperature is changed by varying Nt 
(from 4 to 64). This approach separates temperature and lattice spacing effects and keeps 
the underlying physics scales (i.e. the lattice spacings in temporal and spatial directions, 
respectively at and a,) fixed. 

The calculation of EOS of the quark-gluon system involves derivatives of the bare pa- 
rameters with respect to the physical anisotropy ^ = a.^/at and the spatial lattice spacing 
as- These are the already mentioned above Karsch coeflicicnts, which can be obtained non- 
perturbatively as a by-product of the zero-temperature calculations needed to choose the 
bare parameters. In our scheme, once calculated, the Karsch coefficients can be used for 
all temperatures since they depend only on the intrinsic lattice parameters and not on A^^. 
This allows a straightforward determination of the temperature dependence of the energy 
and pressure, again at fixed lattice spacing. With two or more slightly different values for 
at, a high-resolution sampling of temperatures can be investigated. 

This paper is organized as follows. In section II, we define the anisotropic staggered action 
and derive analytic expressions for the energy and pressure for the studied system. In section 
III we describe the scale-setting techniques. Section IV investigates the phase transition for 
the staggered fermions as we change the temperature in small steps. In section V we present 
the technique used to calculate the Karsch coefficients and the numerical results for them. 
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Section VI contains the numerical results for the EOS for two anisotropies, ^ = 4.0(1) and 
^ = 4.8(3). In Section VII we examine our data for evidence for improvement of the flavor 
symmetry, when due to the anisotropy at becomes sufficiently small. Our conclusions are 
given in section VIII. 



II. ACTION AND GENERAL ANALYTIC FRAMEWORK 



In this section we define the action we use and derive the analytical form of the EOS 
using the derivative method. 

We work with an asymmetric lattice in Euclidean space with notation for the spatial 
lattice spacing and temporal lattice spacing a^. Our calculations are based on the QCD 
action S — Sg + Sp, where the gauge part is the standard anisotropic Wilson action [9]: 



Sg^- 



Co 



x,s>s 



(1) 



and the fermion part is the standard staggered action [10] with anisotropy introduced in the 
spirit of [9] : 



(2) 



In the above definitions Pgs' and Pgt are the space-space and space-time plaquette vari- 
ables. The bare anisotropy parameter = ^ at the tree level, and Jps are the temporal 
and spatial parts of the staggered Dirac operator, i/f and are the bare speed of light 
parameters and Nc = 3 is the number of colors. 

After integrating out the fermion fields in the path integral explicitly, the fermion action 
effectively becomes: 



Nf 



(3) 



where Nf = 2 is the number of fermion flavors and M is the fermion matrix which has the 
form: 



So 



(4) 



We choose Vs = 1 at the expense of rescaling the quark flelds and in an actual simulation ut 
is tuned so that the relativistic properties of the action are restored. 

To determine the energy density e{T) and pressure p{T) as functions of the temperature 
T we use the thermodynamic identities: 



e{T) = - 
p{T) = T 



1 ainZ 



Vs d{llT) 
dlnZ 



V, 



(5) 
(6) 
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where the partition function is Z — J dt/jdt/jdU exp(— 5"), the volume is Vs — N^a^ and 
1/T — Ntttt- Ng and Nt are respectively the number of spatial and temporal lattice sites. 
This way we have: 



dS 



a., , dS 



) 



SmNta^.at'da, 



(7) 
(8) 



where the angle brackets denote averaging over the gauge ensemble. 

The physical anisotropy is defined as ^ = a^/a^ It is convenient to choose ^ and to 
be the independent variables in this relation. This allows Eq. (7) and Eq. (8) to be written 
only in terms of derivatives of them via the transformations: 



_d_ 
dat 
_d_ 
da. 



d 



at 



da.. 



+ 



1 ^ 
at 



(9) 
(10) 



Thus the expression for the pressure p{T) becomes: 

dS 



^da 
dS 



SN^Nta^at 

e _ 

3N^Nta% ^ dC 

e{T) as 



, 1 ,dS 



(11) 



)- 



dS 



smNta^at ' das 



dS 



3 SN^Nta^at ' da, 

To simplify the analytic expressions, in the explicit form of the derivatives of the action S 
we use the following normalization notations: 



(ReTr[P,,,]) 
(ReTr[P,,]) 



(E.,.>yReIV[P.y(x)]) 

(S..sReTr[P,,(x)]) 

^N^s^tNc 
(Tr [M-i]) 

(Tr[^,M-i]) 

{Ti[IpsM-^]) 



(12) 
(13) 
(14) 
(15) 
(16) 



The equations for the energy density and pressure at a given temperature, Eq. (7) and 
Eq. (11), are not corrected for the zero temperature divergent contribution, which simply 
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should be subtracted. This subtraction is trivial and from here on the formulae will assume 
that e{T) and p{T) have that correction. Dividing Eq. (7) and Eq. (11) by T'^ (i.e. mul- 
tiplying them by N^af) and using the notations from above we obtain the following final 
formulae: 



£(r) 

rp4 



3N^ 



3NfN 



+ 



+ 



4e 



dlTLj 



] (ReIV[P,t]) 



as' 



(17) 



rp4 



e{T) ^ a^Nt 




+/9 



(ReTr[P,,,]) 




da. 



+ 



4^ 



dm 



da. 



da. 



(ReTr[P,,]) 



— ^ duf 



(18) 



In order to be able to calculate numerically Eq. (17) and Eq. (18) we need to measure 
all the lattice observables in the above equations and to determine the values of the Karsch 



coefficients ^ 



dnif 



dvt 



' da, 



da, 



dnif 
^' daa 



and 

^ das 



III. SIMULATIONS AND SCALE SETTINGS 



For the purpose of our simulations we implement the R algorithm [11] with step-size 
At — 0.005 and stopping condition 10"^. Figure 1 shows that our choice for the step-size 
allows us to measure physical quantities with an error due to finite step-size smaller than 
2%, and that we are running in the stable regime of the R algorithm. For the spectrum 
measurements we use box sources of size 2 and local sinks for all runs. 

Our simulations examine the phase transition and the thermodynamic properties of QCD 
for volumes, temperatures, spatial lattice scales and quark masses similar to the already used 
in the Nt — 4, 2-flavor thermodynamic studies on isotropic lattices [5]. We adjusted the 
bare parameters in the action so that the resulting physical anisotropy ^ 4, while ~ 0.3 
fm and m^/mp ~ 0.3, which allowed the critical Nt ~ 16. 
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Another important step in tuning the bare parameters is choosing ut such that the 
relativistic properties of the anisotropic staggered action are restored. The velocity of hght 
Cts is defined through the meson dispersion relation: 



El 




I ^2 p2 J_ 
^ Ws-' s,lat 2 ' 



E 



where E and P are the energy and the momentum of the meson, subscripts "lat" and "phys" 
refer to a quantity in lattice or physical units, and s and t whether it is measured in the 
spatial or temporal direction. We tune z/j, so that the velocity of light Cfs(-Ps,phys) ~ 1- The 
velocity of light is calculated for the tt propagating in the temporal direction with a non-zero 
momentum for three valence values of this parameter, z/™' = 0.8, 1.0 and 1.2 (the dynamical 
value is vf^^ — 1.0). Figure 2 demonstrates that the choice of = 1.0 gives a velocity of 
light closest to 1.0 for the set of bare parameters (3 = 5.3, = 3.0 and rrif = 0.008. 

In Table 1 we compare masses measured in the temporal and spatial directions for various 
combinations of z/f-^" and z/™'. This shows that a 20% change in i/f-^'^ has only a small effect 
compared to a similar change in vj^. For the masses we measure the essential contribution 
comes from i/J^'. 

Tabic 11, Table III and Table IV list all the zero temperature scale setting runs and the 
finite temperature runs that we have done. The anisotropy ^ for each of the zero temperature 
runs is calculated from the ratio of the p masses in the spatial and temporal directions. For 
runs 1-5 and 7, Table II, ^ is determined from the matching of the static potentials as well 
(Figure 3 illustrates the matching technique [9]). The comparison between the two methods 
can be done examining Figure 4, which shows that they give reasonably close results. 

The quality of our data for all runs from Table 11 can be judged by studying the effective 
mass plots for n and p. On Figure 5 through Figure 8 we show some typical effective 
mass plots. The effective mass at a given time or space slice is calculated from the values 
of the correlators at 2 neighboring points for tt and 4 for p, in order to determine all the 
parameters in the respective one and two cosh fitting forms. Ideally after some minimal time 
or space slice the effective mass plots should exhibit a plateau. For some effective mass plots 
determined from the spatial correlators the quality of the plateau is not high, especially for 
the runs with very coarse lattice spacing, which means that for those runs there are larger 
errors on the meson masses determined from the fits to all data points. 

All the data which we used in the Karsch coefficients and the EOS determination is given 



The finite temperature runs from Table III and Table IV correspond to two sweeps 
through the phase transition for two different anisotropics ^ = 4.0(1) and ^ = 4.8(3), with 
corresponding a, of 0.34(1) fm, and 0.354(9) fm and (m^/mp)*^"**""""' of 0.33(1) and 0.325(8), 
respectively. 

For each sweep through the transition region the temperature is changed only by changing 
Nt, while all other parameters are kept constant. We want to stress the fact that there are 
no scale changes between the different finite temperature runs from a group with a given 



in Table V. 



IV. PHASE TRANSITION 
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anisotropy and that the scales change minimally between the two groups of runs belonging 
to the two different anisotropies. 

Figure 9 shows the temperature dependence of {ipip) in the critical region. From the data 
we can estimate Tc ~ 150 — 160 MeV. An interesting observations is that the shape of the 
transition is comparable in sharpness with the phase transition obtained in previous isotropic 
calculation [12,13] also shown on the same figure. The isotropic data is from a dynamical 
staggered calculation with two fermion flavors on 16^ x 4 volume and m/ = 0.025. The 
scale used to calculate the temperature in the isotropic case is from [5]. The differences 
between our anisotropic result and the isotropic one we attribute to the scaling violations 
in the latter which we have avoided in our fixed parameter scheme. 



V. KARSCH COEFFICIENTS 



To determine the EOS we need to know the Karsch coefficients which are involved in 
the analytic expressions Eq. (17) and Eq. (18). The values of these derivatives can be 
calculated using the physical quantities that we measure for each zero-temperature run : a^, 
^, Rt = (m^/m2)t^'"P°'-^i and Rst = (m^/m^)«P^«^7(m^/m2)*^'"P°'-^i. We consider the bare 
parameters ^o, P, and Ut to be functions of the above physical quantities, which allows 
those functions to be expanded in Taylor series around the physical quantities of a selected 
zero-temperature run as follows: 



A^oit as: Rt: Rst: {Ci}) = CiA^ + CsAa, + CsARt + c^ARst + ■ 

A(3(^, as: Rt: Rst: {di}) = diA^ + daAa, + d^ARt + d^AR^t + 
Am/(^, as, Rt: Rst: {ei}) = eiA^ + e-iAag + esAi?* + e^ARgt + • 
Avtii, as: Rt: Rst: {fi}) = fi^i + /2Aa, + f^ARt + UARst + • 



(19) 
(20) 
(21) 
(22) 



si 



where A^o = io-i'o: A/? = f^-f^': ^rrij = mf-m'f, Aut = i^t-^v = ^-^': Aa, = a, 
ARt — Rt — R't, ARst — Rst — R'st- definitions the primed quantities refer to the 

selected run around whose physical quantities the Taylor expansion is done. The derivatives 
Ci, di, Ci and fi, i — (1, . . . , 4), are defined as: 



C2 C3 C4^ 

di d2 ds d^ 

ei 62 63 64 

V/i /2 /a fij 
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as 


a. 


,Rt,Rst 


das 


^,Rt,Rst 


aso 

dRt 


£,,as,Rst 


d^o 
dRst 




,Rt 




a/3 


a. 


jRtfRat 


d§_ 

das 


i,Rt,Rst 


df3 
dRt 


£,,as,Rst 


df3 
dRst 


S,as 


,Rt 




dmf 






dmf_ 




dmf 




dmf 








9^ 


a. 


,Rt,Rst 


das 


£.,Rt,Rst 


aRT 


^,as,Rst 


dRst 


S,as 


,Rt 


\ 


dut 

as 


1. 


,Rt,Rst 


dut 
das 


£.,Rt,Rst 


dut 
dRt 


i,as,Rst 


dut 
dRst 




,Rt 1 



(23) 



The Karsch coefficients, which are involved in the EOS, are the first two columns of the 
matrix of derivatives above. We assume that we can make linear fits to Eq. (19) through 
Eq. (22) for each zero-temperature run and minimize the %?, i = (1, ... ,4), for all of the 
zero-temperature runs at the same time. The xf 's for the four fits are: 
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x?({ci}) = E[AC - Ac(r, <, Ri, r:„ {c,}ry^{Ao 

r 

xim}) = - ^p^ic, <, ri, r:„ {d,})ry,{Ap^) 

r 

r 

xim) = - A^[(r, a:, ri r:„ {mVcr'^iAi^D. 



(24) 
(25) 
(26) 
(27) 



In the above expressions the sums are over r, which labels each zero-temperature run and 
the bare parameters and physical quantities associated with it. This labeling is not the same 
as the numbering of the runs in Table II, where each run number refers to a specific set of 
dynamical bare parameters. Here the subscript r labels a specific set of bare parameters 
which instead of uf^"^ has the valence value of that parameter, since as we already showed, 
the uj^^ has the dominant contribution to the measured physical quantities. Hence in all 
formulae in this section the notation stands for the valence value of that parameter. 

All expansions are taken around a given selected run, whose label is not shown. The 
minimization of xl for the first fit, Eq. (19) for example, leads to a matrix equation of the 
form 



AC^V, 



(28) 



where 



E 



ARIAC 



E,. 



A{'-A< 




AS,^ARl 




-^^CA^S) 






^?(AC) 


AoJAoJ 




AaJARJ- 




AalARst 




^?(A?S) 


<^?(A«S) 


t2(a?;) 




ARlAal 


Er 


AR^ARl 


Er 


AR^ARl, 






^?{Ae) 


^?(A«S) 




ARLA< 




ARl^ARl 








^^(A€S) 


<^?(A€S) 


<^?(A€;) 


1 



V ^ 



aliAil) 



\C4/ 



which we solve for C. In a similar way we find matrix equations for the rest of the derivatives 
fi-om Eq. (20), Eq. (21) and Eq. (22). 
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To apply statistical analysis on our data from the zero-tempcratTirc runs, we divide 
the data into a set of jackknifc blocks. The numerical procedure for the minimization of 
the Xi functions can not be applied straightforwardly for the equations Eq. (24) through 
Eq. (27) since the standard deviations of the bare parameters, (T^(. . . ), are not known from 
the beginning. Instead we employ an iterative scheme which consists of the following steps: 

• Start by guessing initial values for all af{. . . )'s. 

1. Determine the Karsch coefficients by solving the matrix equation: Eq. (28), and 
the similar equations derived from minimizing Eq. (25) through Eq. (27). 

2. Using the values of the Karsch coefficients from step 1 calculate numerically the 
linear part of the functions Eq. (19) through Eq. (22) on each jackknife block of 
data and by statistically analyzing them find new values for all (T^(. . . )'s. 

• Repeat steps 1 and 2 until the numerical result for the Karsch coefficients converges. 

The success of this scheme depends on how well the functions Eq. (19) through Eq. (22) 

can be approximated by the linear part of the Taylor expansion, which is a measure of how 
"close" the physical quantities measured from each run are to the quantities of the selected 
run around which the expansion is made. To trust the consistency of the iterative scheme 
we checked the two step procedure with variety of random initial guesses for the aj.{. . . )'s, 
which reproduced the same final results. 

The numerical results for the Karsch coefficients from expansion around runs 7 and 8 
from Table II, obtained via the the method described above, are summarized in Tables VII 
and VIII. The quoted errors are calculated using the jackknife method. Our results show 
larger errors on the Karsch coefficients which are derivatives with respect to than the 
errors on those coefficients that are derivatives with respect to ^. A possible explanation 
of that difference could be that the specific parameter space that we explored in our zero- 
temperature runs does not allow a better resolution of some of the Karsch coefficients either 
because it is too limited (we need more runs and more statistics on each of them to improve 
the quality of mass fits) or because the "points" in that space (the zero-temperature runs) 
are not distributed in a favorable way around the run around which we are making the 
expansion, or both. 

VI. EQUATION OF STATE 

In the previous section we described the procedure which allows us to calculate the 
Karsch coefficients needed to determine the QCD equation of state (Eq. (17) and Eq. (18)). 

As stressed before in Section IV we have two groups of finite temperature runs listed in 
Tables III and IV , for each of which we are changing the temperature by only varying and 
keeping the underlying physics scales fixed. Figures 10 and 11 show the numerical results 
for the energy density and pressure for both groups of runs corresponding to anisotropies 
^ = 4.0(1) and 4.8(3). The data is normalized to the continuum Stefan- Boltzmann values 
of the EOS for an ideal relativistic gas for SU (Nf.) color with Nf quark ffavors, which are 
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|| = ^(^c-l + 4^ciV/)^ 12-1725 

and 

PSB £SB 



The errors on the pressure are significantly larger than the errors on the energy due to the 
large errors on those of the Karsch coefficients which are derivatives with respect to and 
are involved only in Eq. (18) for the pressure. The comparison with the free lattice theory 
(squares) gives an explanation of the prominent drop off of e and p in the high temperature 
sector - simply a consequence of the lattice high momentum cut-off. The high momentum 
mode contribution to the EOS becomes dominant with the increase of the temperature, 
which means that at a coarse lattice spacing a high proportion of relevant modes are simply 
excised. Including improvements to the spatial parts of the staggered fermion action would 
be a natural step to reduce those lattice artifacts for high temperatures. 

Our results for the EOS are comparable with the isotropic case [14] in the temperature 
region up to 0.3 GeV for which the cited reference has data. The errors on the energy density 
are comparable with the errors in that isotropic study, or smaller with enough statistics, on 
the other hand the errors on pressure are larger due to the reasons stated above. 



VII. NOTE ON FLAVOR SYMMETRY IMPROVEMENT 

In the continuum and chiral limits the spontaneous symmetry breaking of S'[/(4)^ ® 
SU (4)y in the staggered action, yields 15 Goldstone pions. On the lattice the violation of the 
flavor symmetry leaves us with the remnant U {1) a®U {l)v and only one true Goldstone pion. 
The local pions in the staggered formulation, which fall into 7 irreducible representations, 
are not degenerate any more due to the O(a^) flavor symmetry breaking. However the 
introduction of anisotropy on the lattice makes the lattice spacing in the temporal direction 
much smaller than the spatial one and hence we expect to see an improvement in the flavor 
symmetry. 

We choose A^^ = {ttIt^^ — 'mTr)/mp, where 112 is the second local staggered pion, as a 
quantitative measure of the flavor symmetry breaking in the spatial and temporal directions. 
The data in Table VI shows that in the temporal direction for all runs is smaller than 
its value in the spatial direction, which means that we are seeing improvement of the flavor 
symmetry as at becomes finer. Especially for run 4, the n and 712 look virtually degenerate. 

We expect that the anisotropy has a similar effect on the rest of the pions, although we 
have not investigated numerically how the various mass splittings between them are affected 
by the decrease of a*. 



VIII. CONCLUSIONS 

We have studied the thermodynamic properties of full QCD with 2-flavors of staggered 
fermions on anisotropic lattices. In our calculations we have employed a fixed parameter 
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scheme in which we keep the bare parameters constant and change the temperature by 
varying only the number of the temporal slices Nf. This allowed us to study the phase 
transition for staggered fermions with fixed physics scales. It appears to be comparably as 
sharp as the transition in the isotropic case. 

We have calculated non-perturbatively the Karsch coefficients from series of zero- 
temperatures runs and applied them in the determination of the EOS. Those of the Karsch 
coefficients which are derivatives with respect to have significant errors which are most 
probably due to the limited set of data used in their determination. They respectively give 
rise to large uncertainties in the calculation of the pressure. 

The high temperature behavior of the quark-gluonic system was found to be strongly 
influenced by the underlying lattice cut-off, which gives a maximum temperature at which 
our anisotropic EOS should represent continuum physics. However the fixed parameter 
scheme combined with a spatially improved anisotropic staggered action might give a much 
better result in the high temperature region. 

The anisotropic approach naturally reduces the finite lattice spacing errors associated 
with at and accounts for an improvement of the flavor symmetry for particles propagating 
in the temporal direction. 

It is interesting to mention that our results do not show a pronounced negative pressure 
problem in the confined phase as it has been found in previous EOS calculations using the 
derivative method with perturbatively calculated Karsch coefficients. However, considering 
the generally large statistical errors on the pressure in our calculation, we could not entirely 
exclude the possibility of such a problem being unveiled at low temperatures in a calculation 
with reduced statistical errors. 
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TABLES 



TABLE I. Meson masses for run 1 and 2, Table II, measured with different valence and dy- 
namical uts. The comparison shows that the main contribution to the masses comes from v^^^. 
Notations S and T stand for spatial and temporal directions of measurement. 



run 


dyn val 


mTT, T 




nip, T 


nip, S 


1 


1.0 1.0 


0.31309(28) 


0.57846(83) 


0.6853(63) 


1.218(29) 


2 


1.2 1.0 


0.31023(42) 


0.57677(46) 


0.6897(54) 


1.257(15) 


2 


1.2 1.2 


0.26761(35) 


0.58299(49) 


0.6030(43) 


1.232(32) 



TABLE II. Parameters of zero-temperature calculations. All runs except run 2 have = 1.0. 
Run 2 has i^f^ = 1.2. 



run 


volume 


traj. 


P 


e 


m/ 


1 


16^ X 32 


5800 


5.425 


1.5 


0.025 


2 


162 X 24 X 32 


5100 


5.425 


1.5 


0.025 


3 


16^ X 24 X 64 


1300 


5.695 


2.5 


0.025 


4 


16^ X 24 X 64 


1400 


5.725 


3.44 


0.025 


5 


16^ X 24 X 64 


3400 


5.6 


3.75 


0.025 


6 


16^ X 24 X 64 


4300 


5.286 


3.427 


0.00394 


7 


162 X 24 X 64 


3200 


5.3 


3.0 


0.008 


8 


16^ X 24 X 64 


3000 


5.29 


3.4 


0.0065 


TABLE III. Parameters of finite temperature calculations with £ 


= 4.0(1). 


All runs have 


= 1.0. 










run 


volume 


traj. 


(3 


Co 


mf 


1 


16^ X 24 


8000 


5.3 


3.0 


0.008 


2 


16^ X 20 


9800 


5.3 


3.0 


0.008 


3 


16^ X 16 


21600 


5.3 


3.0 


0.008 


4 


16^ X 12 


9100 


5.3 


3.0 


0.008 


5 


16=^ X 8 


5500 


5.3 


3.0 


0.008 


6 


16^ X 4 


25900 


5.3 


3.0 


0.008 
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TABLE IV. Parameters of finite temperature calculations with ^ = 4.8(3). All runs have 
4'''' = 1.0. 



run volume traj. (3 mf 

1 16^ X 24 1400 5^29 3^4 0.0065 

2 16^ X 20 2700 5.29 3.4 0.0065 

3 IG^ X 16 8900 5.29 3.4 0.0065 

4 16^ X 12 6200 5.29 3.4 0.0065 

5 16^ X 8 3600 5.29 3.4 0.0065 



TABLE V. All zero-tempcraturc data used to determine the Karsch coefficients from fits to 
Eq. (19) - Eq. (22). The run index in the first column corresponds to the run number of Table II, 
which lists the run parameters. Notations S and T stand for spatial and temporal directions of 
measurement. Rt and Rgt are defined in Section V. 



run 








mjr, S 






Rt 


Rst 


1 


1.0 


1.778(46) 


0.31309(28) 


0.57846(83) 


0.6853(63) 


1.218(29) 


0.2087(39) 


1.080(56) 


2 


0.8 


1.495(21) 


0.37050(41) 


0.56794(39) 


0.8227(69) 


1.230(14) 


0.2028(33) 


1.051(29) 


2 


1.0 


1.822(24) 


0.31023(42) 


0.57677(46) 


0.6897(54) 


1.257(15) 


0.2023(32) 


1.041(28) 


2 


1.2 


2.043(56) 


0.26761(35) 


0.58299(49) 


0.6030(43) 


1.232(32) 


0.1970(28) 


1.137(62) 


3 


0.8 


2.767(41) 


0.2451(10) 


0.6710(50) 


0.3144(34) 


0.870(11) 


0.608(12) 


0.979(27) 


3 


1.0 


3.637(80) 


0.1925(26) 


0.6774(48) 


0.2547(39) 


0.926(13) 


0.571(16) 


0.936(37) 


3 


1.2 


3.273(76) 


0.2094(39) 


0.6732(47) 


0.2716(54) 


0.889(13) 


0.594(12) 


0.966(31) 


4 


0.8 


3.722(41) 


0.2048(11) 


0.7643(24) 


0.2481(28) 


0.9235(58) 


0.681(16) 


1.005(23) 


4 


1.0 


4.295(54) 


0.17811(84) 


0.7699(24) 


0.2210(23) 


0.9493(62) 


0.649(14) 


1.013(27) 


4 


1.2 


4.888(70) 


0.15955(89) 


0.7786(23) 


0.1980(22) 


0.9680(65) 


0.649(14) 


0.997(28) 


5 


0.8 


3.990(34) 


0.20466(65) 


0.8456(27) 


0.2706(18) 


1.0796(74) 


0.5721(80) 


1.072(18) 


5 


1.0 


4.658(35) 


0.17661(76) 


0.8540(20) 


0.2382(13) 


1.1096(88) 


0.5497(80) 


1.078(16) 


5 


1.2 


5.260(46) 


0.15637(69) 


0.8624(22) 


0.21374(99) 


1.1242(90) 


0.5352(74) 


1.099(21) 


6 


1.0 


4.60(21) 


0.07413(25) 


0.36410(77) 


0.2811(39) 


1.294(51) 


0.0695(20) 


1.14(10) 


7 


0.8 


3.57(30) 


0.13273(42) 


0.47488(42) 


0.3777(51) 


1.35(11) 


0.1235(33) 


1.00(17) 


7 


1.0 


4.03(14) 


0.11280(51) 


0.47682(38) 


0.3274(39) 


1.320(49) 


0.1187(29) 


1.099(77) 


7 


1.2 


4.903(75) 


0.09800(44) 


0.47950(36) 


0.2857(34) 


1.401(15) 


0.1177(28) 


0.996(30) 


8 


1.0 


4.80(28) 


0.09517(39) 


0.46121(58) 


0.2821(56) 


1.354(57) 


0.1138(45) 


1.02(10) 
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TABLE VI. The second local staggered pion masses ttt-ttj and = (777,7^2 — m,T^)/mp from 
the zero-temperature runs with parameters given in Table II. The run index in the first column 
corresponds to the run number of Table II. Notations S and T stand for spatial and temporal 
directions of measurement. For large values of the anisotropy ^, in the temporal direction 
is significantly smaller than the corresponding value in the spatial direction, for some runs even 
consistent with zero. 



run 


val 




777,7^2 ) -i 




A T 




1 
1 


l.U 


1 ^VQ/^/lftA 
i. / ( 0(^4D ) 


U.4DUJ3(^o / ) 


i.iuy(^zo j 


U.ziol(^4o J 


U.4oD(_Z4 ) 





U.o 


i.4yo(^zi j 


U.OoUo(^ /Oj 


i.z4o(yz j 


U.ziDOi(^oyj 


u.o4y( /y j 





1 n 
l.U 


i.ozz(^z4 j 


n A TOO / /I \ 
U.4/zy(o4 j 


1 1 V/' 1 oN 

i.i / (iz j 


U.zaoy( IZ>) 


u.4/o(^yo ) 




z 


1 
i.z 


Z.U4o(5Dj 


u.oyyo^oij 


i.iiy(ozj 


U.zioo(54j 


U.4o5(z i ) 


Q 




U.o 


z. /D/(^4i j 


U.z4o4(io j 


U./4dz(/0J 


U.UiUb^OD j 


U.UoD4(o/ j 


Q 



l.U 


o.Do / (^oU j 




u. / ( y(^ii j 


u.uuynDo ) 


nil 1 n\ 
U.iiU(^iU j 


3 


1.2 


3.273(76) 


0.2115(40) 


0.7592(80) 


0.008(11) 


0.0967(89) 


4 


0.8 


3.722(41) 


0.2071(11) 


0.8379(71) 


0.0094(42) 


0.0797(77) 


4 


1.0 


4.295(54) 


0.17931(97) 


0.8418(66) 


0.0054(47) 


0.0757(70) 


4 


1.2 


4.888(70) 


0.1598(10) 


0.8594(69) 


0.0011(59) 


0.0834(71) 


5 


0.8 


3.990(34) 


0.20853(83) 


0.9962(78) 


0.0143(29) 


0.1395(68) 


5 


1.0 


4.658(35) 


0.17954(77) 


1.0125(79) 


0.0123(31) 


0.1428(70) 


5 


1.2 


5.260(46) 


0.15877(75) 


1.0380(84) 


0.0112(33) 


0.1562(74) 


6 


1.0 


4.60(21) 


0.0956(14) 


0.58(11) 


0.0762(53) 


0.168(85) 


7 


0.8 


3.57(30) 


0.1693(14) 


1.30(36) 


0.0969(42) 


0.61(26) 


7 


1.0 


4.03(14) 


0.1410(11) 


1.24(16) 


0.0861(42) 


0.58(12) 


7 


1.2 


4.903(75) 


0.1213(13) 


1.181(79) 


0.0816(54) 


0.501(56) 


8 


1.0 


4.80(28) 


0.1142(14) 


1.54(21) 


0.0674(54) 


0.80(16) 



TABLE VII. Karsch coefficients from fitting the data to the liner part of the Taylor expansion 
around run 7, Table II, with ^ = 4.0(1). The order of the coefficients is the same as in the matrix 

in Eq. (23). Each row is obtained from the fitting procedure independently from the other rows 
and the x^'s per degree of freedom for each fit is respectively 1.7, 1.0, 1.8 and 0.7. 



0.61(6) 


9.6(6.2) 


2.0(1.1) 


0.3(2.0) 


-0.017(7) 


-1.5(1.1) 


0.5(2) 


-0.03(28) 


-0.0062(4) 


0.18(5) 


0.068(4) 


-0.003(15) 


0.04(5) 


-5.6(4.8) 


-0.7(8) 


1.1(1.4) 
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TABLE VIII. Karsch coefficients from fitting the data to the hner part of the Taylor expansion 
around run 8, Table II, with ^ = 4.8(3). The order of the coefficients is the same as in the matrix 
in Eq. (23). Each row is obtained from the fitting procedure independently from the other rows 
and the x^'s per degree of freedom for each fit is 1.4, 0.7, 2.3 and 0.9. 



0.59(6) 


9.2(6.2) 


2.2(1.0) 


0.7(1.8) 


-0.015(6) 


-1.0(4) 


0.58(8) 


0.02(8) 


-0.0050(8) 


0.11(7) 


0.05(1) 


0.02(2) 


0.06(3) 


-5.0(4.1) 


-0.8(7) 


0.8(1.2) 
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FIGURES 



0.155 




0.41 



5e-05 



5e-05 
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0.00025 




0.0001 0.00015 0.0002 0.00025 

(At)' 



FIG. 1. Dependence of physical quantities on the step-size At for volume 8^ x 32, /3 = 5.35, 
= 3.5, vf^'^ = z/™' = 1.0 and nif = 0.006. Trajectories per point about 900. We simulate at 
At = 0.005. 
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FIG. 2. Tuning of the velocity of light cts using the dispersion relation in section III for vr. 



The run parameters are j3 = 5.3, = 3.0, ruf = 0.008, vf^^ = 1.0; Measurements are done at 
f™' = 0.8, 1.0 and 1.2. Square of the spatial momentum in lattice units agPs^phys is plotted on 
the horizontal axis. The choice of z/™' = 1.0 gives velocity of light closest to unity for that set of 
parameters. 
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FIG. 3. Potential matching technique apphed to static potentials measured from Wilson loops 
in the {z, x) and the (z, t) planes for run 5, Table II. On this plot they are shown after they are 
made to match by dividing the abscissa for the potential measured from the (z, x) loops by the 
anisotropy ^ ^ 5, so that 14(^ajn) = Vz{asn). 
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FIG. 4. Scatter plots for zero-temperature runs 1-5 and 7, Table II. The renormalized 
anisotropy ^ is calculated both from the p mass ratio in the spatial and temporal directions and 
from static potential matching. The two methods give the same anisotropy but the errors for the 
potential matching results are larger. 
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FIG. 6. Effective mass plot for vr propagating in spatial direction. Run parameters are as in 
Table II and data is measured at the dynamical value of Vf 
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FIG. 7. Effective mass plot for p propagating in temporal direction. Run parameters are as in 
Table II and data is measured at the dynamical value of Vf 
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FIG. 9. The temperature dependence of (V'^) in the region of Tc. Points from anisotropic runs 
with a common symbol have the same anisotropy and physics scales. The isotropic data is shown 
for comparison. From the critical region we estimate Tc ~ 150 - 160 MeV. 
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FIG. 10. Energy in units of the Stefan-Boltzmann limit and a comparison with the free lattice 
theory (squares). The Stefan-Boltzmann law for a relativistic ideal gas for SU{Nc) color with 
Nf quark flavors in the continuum is ^ = ^{N^ - 1 + jNcNf) ^ 12.1725 for Nc = 3 and 
Nf = 2. Points with a common symbol share the same anisotropy and the same physics scales at 
all temperatures. 
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FIG. 11. Pressure in units of the Stefan-Boltzmann limit and a comparison with the free lattice 
theory (squares). The Stefan-Boltzmann law for a relativistic ideal gas for SU{3) color with Nj = 2 
quark flavors in the continuum gives = fff • Points with a common symbol share the same 
anisotropy and the same physics scales at all temperatures. 
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